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g_J ' Abstract. We search for CP violation in a sample of 4.7 x 10^ Cabibbo suppressed D^ — >■ 

K* K~ 7r+ tt" decays. We use 470 fb"^ of data recorded by the BABm. detector at the PEP-II 
asymmetric-energy e^e" storage rings running at center-of-mass energies near 10.6 GeV. CP 

pH ' violation is searched for in the difference between the T-odd asymmetries, obtained using triple 

product correlations, measured for D" and D decays. The measured At violation parameter 
is At = (1.0 ± 5.1stat ± 4.4syst) x 10"^ 
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In Standard Model, CP violation arises from Kobayashi-Maskawa phase in Cabibbo- 
Kobayashi-Maskawa quark mixing matrix|]4i2j. Theoretical attempts to predict the effect of CP 
_ violation in Cabibbo suppressed charmed decays have been made in the past[3j, obtaining a limit 

Q ■ of 0.1% not excluding even 1% effects. The same paper suggests that this limit can be lowered 

f^ I by at least one order of magnitude by oscillations, which have been recently observed [H |5]. 

CP violation in charm decays can be exploited by many New Physics models [6] both at 
tree and one-loop level; among these the latter expect a CP violation asymmetry at the order 
'k^ [ 0(10"^), which is now the level of experimental sensitivity [7]. 

^ ■ We make use of T-odd correlations [8] to build a T odd observable: assuming CPT theorem, 

. . , CP violation is straightforward once T violation is found. One way to build a T odd observable 

rely on the mixed product vi ■ {v2 x v^), where each Vi is a momentum or a spin. A non-zero 
triple product correlation is then evidenced by the asymmetry 

^ ^ Tjvi ■ {V2 X ^3) > 0) - r(^i • {V2 X ^3) < 0) 

^ ~ r(i;i • {v2 X V3) >0)+ r(7;i • {v2 x v^) < 0) ' 

where T is the decay rate of the process. There is however a technical complication due to strong 
phases, which can fake this signal. The true T violation observable is then 

At = -{At- At), 

where At is the charge conjugate of At, in which by definition the weak phase changes its sign, 
while the strong does not. This observable can be built in the D^ — )■ K^ K~ 7r+ vr" decays 
defining Ct = 'Pk+ ' {Pn+ ^ Ptt- )) using the momenta pi of the final state particles in the D^ rest 
frame, and taking (Ct = Pk- " (Pn- ^ Pn+)) 

_ T{Ct > 0) - TJCt < 0) T _ TJ-Ct > 0) - r{-CT < 0) 

^ ~ r{CT > 0) + r(CT < 0) ^ ~ r(-CT > o) + r{-CT < o) ' 



The reaction [9] 



e e 



X D*+: D*^ 



vr+D"; D^ 



K+K-TT+TT- 



where X indicates any system composed by charged and neutral particles, has been reconstructed 
from the sample of events having at least five charged tracks. We first reconstruct the D^ 
candidate: all K~^K^it^-it^ combinations assembled from well-measured and positively identified 
kaons and pions are constrained to a common vertex. To reconstruct the D*~^ candidate, we 
perform a vertex fit of the D^ candidates with all combinations of charged tracks having a 
laboratory momentum below 0.65 GeV/c (vr^) with the constraint that the new vertex is located 
in the interaction region. We require the D^ to have a center-of-mass momentum greater than 
2.5 GeV/c: this requirement removes D^ coming from B decays. We observe a contamination of 
the signal sample from D^ — )• K^K~Kg, where Kg — )• tt^tt~. The tt+tt" effective mass shows, 
in fact, a distinct Kg mass peak, which can be represented by a Gaussian distribution with 
a = 4.20 lb 0.26 MeV/c^, which accounts for 5.2% of the selected data sample. We veto Kg 
candidates within a window of 2.5 a. This cut, while reducing to negligible level the background 
from D^ — )• i^+ K~ Kg, removes 5.8% of the signal events. 

Defining the mass difference Am = m{K'^ K^ tt^ tt^ irf) — m(K^K^TT^iT~), Figure H^a) 
shows the scatter plot m{K~^K~7r~^7r~) vs. Am for all the events. Figure W(Jo) shows the 
m.(K~^K~ Tr~^Tr~) projection, Fig. [^c) shows the Am projection. 
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Figure 1. (a) m{K^K^TT^TT^) vs. Am, for the total data sample, (b) m{K^K^'iT^TT^) 
and (c) Am projections with curves from the fit results. Shaded areas indicate the different 
contributions. The fit residuals, represented by the pulls, are also shown under each distribution. 



We perform a fit to the m.{K~^K~TT~^TT~) and Am distributions, using a polynomial 
background and a single Gaussian. The fit gives cr^jo = 3.94 zt 0.05MeV/c^ for the D^ and 
CAm = 244 ± 20keV/c^ for the Am. We define the signal region within ±2a£)0 and ib3.5(TAm- 
The total yield of tagged D^ mesons in the signal region is approximately 4.7 x 10^ events. 

The D^ yields to be used in the calculation of the T asymmetry are determined using 
a binned, extended maximum-likelihood fit to the 2-D {m{K~^K~7r^TT~), Am) distribution 
obtained with the two observables m.{K^K^iT^Tr^) and Am, in the mass regions defined in 
the ranges 1.825 < m{K+K--K+TT-) < 1.915 GeV/c^ and 0.1395 < Am < 0.1545 GeV/c^ 
respectively. Events having more than one slow pion candidate in this mass region are removed 
(1.8 % of the final sample). The final 2-D distribution contains approximately 1.5 x 10^ events 
and is divided into a 100 x 100 grid. 



The 2-D {m{K^K vr+vr ), Am) distribution is described by five components: 

(i) True D^ signal originating from a D*^ decay. This component has characteristic peaks in 
both observables m^K^K^ir^iT^) and Am. 

(ii) Random vr^ events where a true D^ is associated to an incorrect tt^, called D^ peaking. 
This contribution has the same shape in m{K~^K~TT~^TT~) as signal events, but does not 
peak in Am. 

(iii) Misreconstructed D^ decays where one or more of the D^ decay products are either not 
reconstructed or reconstructed with the wrong particle hypothesis, called Am peaking. 
Some of these events show a peak in Am,, but not in m{K~^K~Tr~^TT~). 

(iv) Combinatorial background where the K~^ , K~ , vr^, vr^ candidates are not fragments of 
the same D^ decay, called combinatoric. This contribution does not exhibit any peaking 
structure in m,{K~^ K~ tt^ tt~) or Am. 

(v) D^ — )• K~^ K~ vr"*" 7r~ vr"*" contamination, called Df. This background has been studied on 
Monte Carlo (MC) simulations and shows a characteristic linear narrow shape in the 2-D 
(m(K+K~7r+7r~), Am) distribution, too small to be directly visible in Fig. [TJa). 

The functional forms of the probability density functions (PDFs) for the signal and 
background components are based on studies of MC samples. However, all parameters related 
to these functions are determined from two-dimensional likelihood fits to data over the full 
m,{K~^K~TT~^Tr~) vs. Am region. We make use of combinations of Gaussian and Johnson SU [11] 
lineshapes for peaking distributions, and we use polynomials and threshold functions for the non- 
peaking backgrounds. The results of the fit are shown in Fig. [TJ The fit residuals shown under 
each distribution are represented by Pull = [N^ata — ^fit)/V^data- 

According to the D*~^ tag and the Ct variable, we divide the total data sample into four 
subsamples, defined in Table [TJ These four data samples are fit simultaneously to the same 
model. The signal event yields are given in Table [TJ 

Table 1. Definition of the four subsamples and the event yields from the fit. 



Subsample 


Events 


(a) D'>, Ct>0 


10974 ± 117 


(b) DO, Ct <0 


12587 ± 125 


(c) D', Ct > 


10749 ± 116 


(d) D^, Ct<0 


12380 ± 124 



We validate the method using e'^e~ — )• cc MC simulations, where D^ decays through the 
intermediate resonances with the branching fractions reported in the PDG [13j. We obtain a T 
asymmetry At = (2.3 it 3.3) x 10~^, consistent with the generated value of 1.0 x 10^'^. 

To test the effect of possible asymmetries generated by the detector, we use signal MC in 
which the D^ decays uniformly over phase space. In this case possible asymmetries are generated 
only by the detector efficiency: At = —(1-1 i 1.1) x 10"'^, again consistent with zero. 

To avoid potential bias, all event selection criteria are determined before evaluating At- 
Systematic uncertainties are obtained directly from the data. In these studies the true At and 
At central values are masked by adding unknown random offsets. Removing the offsets: 

At = (-68.5 ± 7.3stat ± 5.8,y,t) x 10"' At = (-70.5 ± 7.3stat ± 3.9syst) x W~\ (1) 

We observe non-zero values of At and At indicating that final state interaction effects are 
significant in this D^ decay. No effect is found, on the other hand, in the analysis of MC samples. 



The result for the CP violation parameter, At, is 

At = (1.0 ± 5.1stat ± 4.4syst) x 10"'. (2) 

The sources of systematic uncertainties considered in this analysis and the estimates of their 
values are derived as follows: 

(i) The PDFs used to describe the signal are modified, replacing the Johnson SU function by 
a Crystal Ball function [12j, obtaining fits of similar quality {asys = 0.2 x 10"'^). 

(ii) As the same as (i), for the peaking background {agys = 0.5 x 10^^). 

(iii) We increase the number of bins of the 2-D (m^K^K^ir^Tr^), Am) distribution to a 
(120 X 120) grid and decrease to a grid of (80 x 80) {asys = 0.2 x 10^^). 

(iv) The particle identification algorithms used to identify kaons and pions are modified to more 
stringent conditions in different combinations {asys = 3.5 x 10~^). 

(v) The p*{D^) cut is increased to 2.6GeV/c and 2.7GeV/c {asys = 1.7 x 10"^). 

(vi) We study possible intrinsic asymmetries due to the interference between the electromagnetic 
e^e~ — )• 7* — 7- cc and weak neutral current e^e~ — )■ Z*^ — t- cc amplitudes. This interference 
produces a D^ /D^ production asymmetry that varies linearly with the quark production 
angle with respect to the e~ direction. We constrain the possible systematics by measuring 
At in three regions of the center-of-mass D^ production angle 6*: forward (0.3 < cos{6*)t)o), 
central (-0.3 < cos(6'*)£,o < 0.3), and backward (cos(6l*)£,o < -0.3) {asys = 0.9 x 10^^). 

(vii) Fit bias: we use MC simulations to compute the difference between the generated and 

reconstructed At {asys = 1.4 x 10"'^). 
(viii) Mistag: there are a few ambiguous cases with more than one D* in the event. We use MC 

simulations where these events are included or excluded from the analysis. This effect has 

a negligible contribution to the systematic uncertainty. 

(ix) Detector asymmetry: we use the value obtained from the MC simulation where D^ decays 
unifromly over the phase space {asys = 1.1 x 10~^). 

In the evaluation of the systematic uncertainties, we keep, for a given category, the largest 
deviation from the reference value and assume symmetric uncertainties. Thus, most systematic 
uncertainties are statistical in nature, and are conservatively estimated. 

In conclusion, we search for CP violation using T-odd correlations in a high statistics sample of 
Cabibbo suppressed D^ — )• K'^K^tt^tt^ decays. We obtain a T-violating asymmetry consistent 
with zero with a sensitivity of ~ 0.5 %. These results constrain the possible effects of New 
Physics in this observable. 
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